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We investigate a scheme that makes a quantum non-demolition (QND) measurement of the exci- 
tation level of a mesoscopic mechanical oscillator by utilizing the anharmonic coupling between two 
beam bending modes. The non-linear coupling between the two modes shifts the resonant frequency 
of the readout oscillator in proportion to the excitation level of the system oscillator. This frequency 
shift may be detected as a phase shift of the readout oscillation when driven on resonance. We de- 
rive an equation for the reduced density matrix of the system oscillator, and use this to study the 
conditions under which discrete jumps in the excitation level occur. The appearance of jumps in the 
actual quantity measured is also studied using the method of quantum trajectories. We consider 
the feasibility of the scheme for experimentally accessible parameters. 

PACS numbers: 63.20.Ry, 03.65.Yz 

I. INTRODUCTION 

Quantum mechanics tells us that the energy of an oscillator is quantized. However, an observation of quantum 
limited mechanical motion in macroscopic objects has not been possible because the energy associated with individual 
phonons is typically much smaller than the thermal energyii^. 

Advances in nanotechnology have enabled experimenters to build ever smaller mechanical oscillators with high 
resonance frequencies and quality factors'^. As an individual phonon energy becomes comparable to or greater than 
fceT, quantum effects begin to appear and it should be possible to realize various quantum phenomena. 

In this paper, we investigate the possibility of observing transitions amongst the Fock states of a mesoscopic 
mechanical oscillator. To do this requires the coupling of the system oscillator to a measurement device that sensitively 
detects the phonon number of the system oscillator but does not itself change the excitation level of the oscillator. 
In the quantum regime it becomes very important to model the precise way that a quantum system interacts with 
any measuring apparatus, as well as with the environment. Specifically, it is necessary to take into account the 
measurement backaction and to design the system-readout interaction so as to allow the best possible measurement 
of the desired observable. We will show that it is possible in principle to take advantage of the non-linear interaction 
between modes of oscillation of an elastic beam or beams to track the state of the oscillator as it jumps between 
number states due to its coupling to the surrounding thermal environment. 

The laws of quantum mechanics tell us that, even in the absence of instrumental or thermal noise, a measurement 
will tend to disturb the state of the measured system. The interaction between the system and the measurement 
apparatus means that while information about the measured observable may be read out from the state of the meter 
after interacting with the system, the quantum mechanical uncertainty in the initial state of the meter leads to 
random changes in the conjugate observable of the system. This backaction noise on the conjugate observable is 
an inevitable result of the very interaction that allows the measurement to take place. It has long been recognized 
that such backaction noise places a fundamental limit on the sensitivity of physical measurementSi^. However, the 
class of measurements known as quantum non-demolition (QND) measurements partially circumvents this problem by 
guaranteeing that the backaction noise does not affect the results of future measurements of the same quantity. The 
idea of a QND measurement is widely discussed in the literature (for example, see^f£i, and^). In a QND measurement 
the interaction Hamiltonian between system and meter commutes with the internal Hamiltonian of the system: an ideal 
QND measurement is repeatable since the backaction noise does not affect the dynamics of the measured observable. 
In this paper, we are interested in a QND measurement of phonon number. The conjugate observable of number is 
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phase, thus, the measurement backaction in our case will result in diffusion of the phase of the mechanical oscillations. 
However, the scheme allows the complete determination of the oscillator excitation level and thus projects the system 
onto a number state in an idealized limit. 

The scheme for the QND measurement of phonon number that we consider uses two anharmonically coupled modes 
of oscillation of a mesoscopic elastic structure. The resonant frequencies of these two modes are different. The higher 
frequency mode is the system to be measured, while the lower frequency oscillator serves as the meter (we refer this 
oscillator as ancilla). The key idea of the scheme is that, from the point of view of the readout oscillator, the interaction 
with the system constitutes a shift in resonance frequency that is proportional to the time-averaged excitation of the 
system oscillator. This frequency shift may be detected as a change in the phase of the ancilla oscillations when driven 
on resonance. We show that this scheme realizes an ideal QND measurement of phonon number in an appropriate 
limit. To measure the phase of the ancilla oscillator we imagine a magnetomotive detection scheme so that the 
actual physical quantity measured is an electric current that couples to the ancilla displacement. Thus our task is to 
understand how the strong measurement of the current, represented by the von Neumann projection scheme on the 
current operator, yields information on the system phonon number, and in turn affects the dynamics of the system via 
the indirect coupling through the ancilla oscillator, and in the presence of the inevitable coupling of the ancilla and 
system oscillators to the environment. This QND measurement scheme where a nonlinear potential provides a phase 
shift to one oscillator that reflects the excitation of the other is analogous to the experiment of Peil and Gabrielse^, 
which demonstrated a QND measurement of the excitation of a single trapped electron. Theoretical discussions of 
such approximate QND measurements of the excitation of an oscillator date back at least to Unruh^ . 

In Sec. ^ we introduce our model and construct the Hamiltonian describing the two oscillators, the magnetomotive 
coupling, and the coupling to the environment represented by baths of harmonic oscillators. For the ancilla oscillator 
displacement to directly indicate the system phonon number the time scale of the ancilla dynamics must be much 
shorter than that of the system. This actually allows us to adiabatically eliminate the ancilla operators to obtain 
dynamical equations for the system alone. Thus, in Sec. IIIII we obtain a reduced master equation for the density 
matrix of the system oscillator, which allows us to focus on the physics of the system dynamics. This adiabatic 
elimination of ancilla degrees of freedom is often considered in quantum optics. However, the adiabatic elimination 
used in quantum optics is at temperature zero, and we need to reformulate the method for finite temperatures. 

Once we know the system dynamics, we next focus on obtaining the experimental outcome. Quantum mechanics 
allows us to determine the state of the system conditioned on the measured current I{t). The von Neumann pro- 
jection postulate says that after a measurement a quantum system in some possibly mixed initial state is projected 
onto the eigenstate corresponding to the measurement outcome. The continual measurement and projection of the 
current I{t) provides accumulating information on the system phonon number, and correspondingly a projection onto 
phonon number states. The theory of quantum trajectorieaSiiSiiiii^ has been developed to deal with such continuous 
measurements. In Sees. IIVI we discuss such quantum trajectory equations for our system. The method leads to a 
stochastic master equation for the system density matrix, where the stochastic component comes from the particular 
value of the measured current at each time, which itself is a stochastic variable since it is the outcome of a quantum 
measurement. These equations of motion for the system conditioned on a particular sequence of measurement results 
allow us to investigate the possibility of tracking the evolution of the system as it jumps between number states due 
to its interaction with the thermal bath. Some details of the formulation of the operators describing the measurement 
current that are needed to derive the stochastic master equation are given in Appendix 1X1 

The main discussion of the physical implications of the model is in Sec. Ivl where we consider the parameters that 
are necessary to observe the oscillator jump between number states. A reader not interested in the details of the 
derivation could read Sec. and following sections after the description of the model in Sec. ^1 Finally, in Sec. IVII 
we conclude with a discussion of the feasibility of the scheme based on current technology and future enhancements. 

II. CONSTRUCTING THE HAMILTONIAN 
A. The model 

In this section we introduce the model system and show how the coupling between the system and ancilla ap- 
proximates a QND coupling in an appropriate limit. We then derive equations of motion that take into account the 
couplings to the environment and the interactions that drive and monitor the oscillations of the ancilla. 

Consider a mesoscopic beam with rectangular cross section. There are two orthogonal flexing modes that are not 
coupled in the linear elasticity theory, but are coupled anharmonically. This coupling exists in nature between the two 
orthogonal flexing modes of single mechanical beam. However, the coupling can also be controlled and engineered: a 
similar coupling of bending modes in two elastic beams has been proposed by YurkeiSi. In this scheme two mesoscopic 
elastic beams with rectangular cross section are connected by a series of mechanical coupling devices. These devices 
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FIG. 1: Schematic of a QND measurement using two coupled mechanical oscillators. 

have the effect of allowing only one type of strain (the longitudinal stretch) to pass to the other beam. In this paper 
we focus on the extent to which such mechanical devices are able to realize a QND measurement and the constraints 
this places on the specifications of the device, and the temperature at which the experiment is performed. 

For convenience, we refer to the system of interest as oscillator and the ancilla as oscillator 1, and the corresponding 
resonant frequencies of the two modes as ujq and uji , respectively. The ancilla is driven at its resonant frequency, and 
a measurement apparatus is attached to the ancilla. The whole structure is kept at a low temperature T such that 
huio ~ ksT, where h — h/2TT is Plank's constant, and is Boltzmann's constant. The oscillators are weakly coupled 
to the environment. Figure ^ shows a schematic of our model. 

B. System Hamiltonian 



Converting the schematic model in the last section into a tractable mathematical model and obtaining the system 
dynamics requires some assumptions and simplifications. 

Firstly, we focus on the anharmonic coupling and the limit in which it satisfies the QND condition. In linear 
elasticity theory, the two flexing modes, which are perpendicular to each other, propagate independently without 
interacting. Beyond the linear approximation these modes are coupled. Expansion of the elastic energy with respect 
to the strain tensor is taken up to second order in the harmonic approximation. By symmetry the coupling between 
the modes first occurs at fourth order, proportional to xf^x\. So we expand the anharmonic terms up to quartic order 
to give 

i/anh = h {Xmxt + ~\iix\^ + hXaixlxl, (1) 

where \ij give the strengths of the nonlinear terms. The first two terms in Eq. are internal anharmonic terms. 
Under the rotating wave approximation (see below), these terms cause a shift in the oscillator resonant frequencies 
and a non-linear phase shift that depends on intensity (a Kerr non-linearity) resulting in rotational shear of the state 
in the phase space of the two oscillators. For the system oscillator, the small shifts in the energy level spacings 
are not important, and can be ignored. The ancilla oscillator is externally driven, and so the nonlinearity of this 
oscillator may become large, for example leading to multistability for large enough drive strengths. We will assume 
that the drive strength is kept smaller than this range, so that again the x\ nonlinearity does not play an essential 
role. However the x'^x\ coupling plays an essential role in coupling the system and ancilla. Therefore, in the interests 
of a straightforward presentation we retain the non-linear coupling given by Aqi and disregard the nonlinearities of the 
system and ancilla internal Hamiltonians. A detailed analysis including non-interacting nonlinearities and detuning 
in Ref. (~) has shown that in the regime of strong damping of the ancilla that we will mostly consider the effect of 
these anharmonic terms will be negligible for small detuning. 

In terms of creation and annihilation operators, the Hamiltonian is now H = Hq + Kinh, with Hq the harmonic part 



Hq = hujoalao + hu!ia\ai, 



and Vanh the anharmonic coupling 
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We have defined the standard raising and lowering operators for the oscillators 



Gi = y/miiVi/2hxi + iy/T/2hm~LJiPi, (4) 

and a| is the Hermitian conjugate of a;. So far we have ignored any coupling of the two oscillators to the environment 
so as to focus on the interaction of the two oscillators. 

In order to perform a QND measurement of alao the Hamiltonian of the oscillators H should satisfy the QND 
condition 

[alao,Ho + V^nh]^0. (5) 

To show that this condition is satisfied in an appropriate limit it is useful to move into an interaction picture with 
respect to Hq. If the frequencies of the two oscillators satisfy — i-^i 3> Aqi and uJi 3> Aqi, then the time-dependent 
terms in t4,nh(t) lead to rapid, small amplitude oscillations of that essentially average to zero over the time 
scales for which the non-linearity Aqi is relevant. If we admit a time coarse-graining over times longer than the 
mechanical oscillation period we may ignore the rapidly oscillating terms, an approximation known as the rotating 
wave approximation (RWA) . Another intuitive explanation for the rotating wave approximation is that the condition 
uq — ui 3> Aoi means that the differences in energy are so large that the energy non-conserving transitions are strongly 
suppressed. 

Disregarding the energy non-conserving terms in the Hamiltonian and then absorbing constant corrections to the 
system and ancilla oscillation frequency into the definition of loq and uji , we obtain 

V^Tit) = hXoialaoaW (6) 

The constant term has been disregarded since it merely provides an overall phase. Note that having made the rotating 
wave approximation the anharmonic coupling term commutes with the observable ajap, and so a QND measurement 
can be achieved under the condition ujq — lui ^ Aqi^^. Returning to the Schrodinger picture the Hamiltonian H now 
can be written as 

jjHWA _ fj^^u^^Q^ -|_ fujj-^^a\ai + S,AoiaJaoa|ai. (7) 

In the above rotating wave Hamiltonian an excitation of the system oscillator leads to a frequency shift of the 
ancilla oscillator. To detect the system excitation level, the ancilla is driven on resonance and the phase shift of 
the oscillations is measured. The driving of the ancilla may be written in terms of term in the Hamiltonian in the 
Schrodinger picture 

^^drivc = 2/l£' (fli + a\] COSUJit, (8) 



where the parameter E is used to characterize the strength of the drive. In the interaction picture using the RWA for 
uji > E, we get 



-"drive ~ 



(ai + al) . (9) 



Now we add the coupling of thermal baths to the system and ancilla. We employ a standard technique and model 
the thermal baths (the surrounding environment) as an infinite number of harmonic oscillators. The thermal baths are 
linearly coupled to the system or ancilla by coordinate-coordinate coupling, i.e., AijXiXj where Xi is the system 
or ancilla coordinate, Xj is the coordinate of an oscillator in the bath, with the index j corresponding to different 
bath oscillators. We will again use the rotating wave approximation for the coupling since the couplings are weak. 

The nature of the coupling with the measurement instrument depends on the measurement scheme. Here we adopt 
a magnetomotive detection scheme suggested by Yurke et a/, i5,i6^i7 _ ^ metallized conducting surface on the ancilla 
oscillator develops an electromotive force across it due to a perpendicular magnetic field and the oscillation of the 
beam. The voltage developed depends on 

V = IB% (10) 

where V is the voltage, B is the magnetic field, the conductor is of length / and xi is the displacement of the beam 
from its equilibrium position. Depending on the resistance R of the conducting strip and remainder of the circuit this 
will result in a current that is then measured. 
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In order to quantize this measuring device we follow the standard practice in quantum electronics and model this 
resistance by a semi-infinite transmission linei^. This model has been used in the context of mechanical measurements 
(seei^iSli) and is in fact mathematically the same as the 'Rubin modeliSiiSS. This is certainly a simplified model of 
an actual detection circuit which essentially assumes that the noise in the circuit is broad-band and Gaussian. The 
transmission line will be considered to be at a temperature corresponding to the effective noise temperature of the 
detection circuit and this noise will affect both the sensitivity and the heating of the ancilla and system oscillators. 
More realistic quantum mechanical models of amplifier circuits can be found in^'^ for example. Our final model for 
the QND set-up will be fairly robust to the precise detection circuit. The resulting current operator is 

/«5]fed.„ + fcL, (11) 

n 

where 6d.n are the lowering operators for the modes of the transmission line and the proportionality constant, which 
is not important for our results, depends on the circuit parameters. For a linearly coupled system-bath measurement 
within the rotating wave approximation, the Hamiltonian describing the coupling between each measurement current 
mode and the ancilla is proportional to 6^ „ai -I- bd,na\- The coupling to the thermal bath modes has the same 
mathematical structure. In the rotating wave approximation, the difference between a coordinate-coordinate coupling, 
and a momentum-coordinate coupling can be absorbed into the definition of the phase of the various raising and 
lowering operators. As is usually done, for later convenience we will include a phase factor of 7r/2 so that the coupling 
to the baths and measurement current takes the latter form. 
The final Hamiltonian for our model is then 

oo 

H = hujQalao + HE (^ai + + X] '^s,n&l.„&s,n 

s n 

+ hXoialaoa\ai + ih (e'^ao - OaJ^ + ih (v^'ai - Ta\^ + ih (^L>^ai - D (t) , (12) 

where T,D,Q have the form, 9s{'^n)bs,n, and the index s denotes the three different baths: the thermal bath 
coupled to the system (0), the thermal bath coupled to the ancilla (1), and the measurement bath coupled to the 
ancilla (d). The strength of the coupling to the bath modes is given by the coefficients gs {ujn)- Later, we will derive 
the relationship between these coefficients and the corresponding oscillator damping rates or quality factors. 

III. DYNAMICS OF THE SYSTEM 

We use the dynamics described by Eq. (|12|l to understand the measurement process. In this section we first find a 
master equation that describes the evolution of the system and ancilla alone without explicitly describing the state of 
the environment. Secondly we further simplify this equation by making use of the difference in time-scales between 
the system and ancilla to obtain a master equation for the system oscillator alone by means of adiabatic elimination. 
This allows us to study the effect of the QND measurement coupling on the system. 

A. Master equation 

We develop a master equation for the density operator of the system alone by integrating out the bath degrees 
of freedom. Because we are interested in high-Q oscillators weakly coupled to the baths we employ a rotating wave 
approximation and the Markov approximation that the memory time scale of the bath is short. In this regime the 
rotating wave master equation accurately describes the dynamics on time-scales longer than an oscillation period, and 
the resulting master equation preserves the positivity of the density matrix. The derivation of such master equations 
is widely discussed in the literature, see for example Carmichael^, Walls and Milburn^ and Caldeira and Leggett^. 
Note that Caldeira and Leggett do not make the rotating wave approximation (which we adopt from the quantum 
optics literature) but instead make a high-temperature approximation; the two equations agree in the overlap of their 
domain of validity (high temperature and weak coupling to the bath). However, the Caldeira-Leggett master equation 
can only be guaranteed to preserve the positivity of the density operator in the limit of high temperature. 

Assuming that the environment and measurement baths are in thermal equilibrium, the master equation for the 
reduced density operator p describing the state of the system and ancilla in the interaction picture takes what is 
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known as Lindblad form 

dp i 
'dt ^ h 



HE 



(^a\ + ai^ + hXoialaoa\ai, 



+ v{Nq + 1)V [ao] p + vNqV [ajj p + k{Ni + l)V [oi] p + kN{D \a}^ p 



(13) 



where 



xp - 



px'^x, 



V [x] p ~ 2xpx^ 

and A'i are the Bose-Einstein factors at frequency loq for A'o and loi for A^i. The first term in Eq. (|13|l involving 
the commutator describes the coherent driving of the ancilla oscillator and the non-linear coupling between the two 
oscillators in the rotating wave approximation. The remaining terms describe the dissipative interactions with the 
various baths. The coefficient v is 

= T^Qho (wo) |.gbo (^o)l^ , 

where g^o (wq) is the density of states of the bath coupled to the system at frequency Wp- It can be experimentally 
obtained from the quality factor of the system oscillator Qq as v = wq/SQo- The coefficient k is the corresponding 
damping rate of the ancilla, with contributions 77 from the coupling to the environment and p from the measurement 
bath. Both these rates can be expressed in terms of the bath density of states at frequency lui in exactly the same 
way as for i^, and k = 7] + p. The terms containing a factor (Ni + 1) describe the spontaneous and stimulated emission 
of phonons into the thermal bath while the ones with Ni correspond to absorption of phonons. 

The master equation Eq. H13|l can in principle be numerically integrated. However we will make some further 
approximations in order to derive a master equation for the system dynamics alone and show that in some limit the 
readout system coupling results in the phase diffusion that is required as the backaction for the QND measurement, 
with no extra noise above this quantum limit. To do this we assume that the ancilla is strongly damped. In this 
limit the ancilla relaxes rapidly to a state consistent with the instantaneous system state. As a result its dynamics 
are slaved to the system oscillator and can in fact be eliminated from the equations of motion. Experimentally this 
is the limit in which the displacement of the ancilla directly reflects the system behavior. This adiabatic elimination 
is described in the following subsection. The final result of this analysis is Eq. H25|l below. 

Note that v, k are the widths of the oscillator resonances, and these should be taken into account when assessing 
the validity of the rotating wave approximation. In the presence of the coupling to the baths, the rotating wave 
approximation is only valid if loq — wi is much greater than the line width of the oscillators, i.e., loq — uji 3> v, k. This 
condition can be understood as not allowing non-energy-conserving transfers of a phonon between the two oscillators. 



B. Adiabatic elimination 



For a strongly damped ancilla (k ^ i') the driven ancilla rapidly relaxes to a state that oscillates with a phase 
determined by the current system phonon number. In the interaction picture this is a displaced thermal state: a state 
with variance of position and momentum equal to those of a thermal state but with non-zero expectation values of 
position and momentum consistent with the driving and damping of the oscillator. It will be useful to transform the 
equations of motion in such a way as to make a perturbative expansion around this steady state. The basic idea is to 
transform the origin of phase space such that the ancilla steady state for the transformed master equation is a thermal 
state. This transformation will essentially remove the driving term in the master equation. This is the approach of 
Wiseman and Milburn^® who study adiabatic elimination in a similar model. While they assume zero temperature 
and therefore end up with a perturbation expansion about the displaced ancilla ground state we must generalize their 
techniques to finite temperature. 

Following Wiseman and Milburn, we use the displacement operator, D{a) — exp aa[ — a*ai with a — —iE/n. 
The transformed system state is p = D{a)pD{a)'^ and we may write the master equation for p as 

p = D{a)pD{a)^ 
= -i|a|^Aoi alaa,p 



iXo 



aJaoa|ai, p 



iXoi 



k{Ni + 1) (^2aipa\ — a\aip — pa\aij + uNi (^2a\pai — aia\p — paia\ 
v{No + 1) {2aopal - ala^p - pa^a^^ + z^A^o {^alpao - paoal - aoajp 



(14) 



7 



In this master equation the excitation of the ancilla oscillations leads to a frequency shift of the system oscillator 
described by the first three terms on the right hand side of this equation. The first term is due to the classical mean 
value of the ancilla oscillator energy and is just a constant shift in the system oscillation frequency. We may move 
to an interaction picture at this shifted frequency, the most convenient interaction picture in which to perform the 
adiabatic elimination. The next two terms describe the effect of the fluctuations in the ancilla excitation. The thermal 
bath coupling terms (the last four groups of terms) are the same as before. 

The adiabatic elimination will hold when the terms proportional to n in Eq. (|14|l dominate in the ancilla dynamics. 
Thus, the adiabatic elimination is valid in a strongly damped regime such that 



(15) 



We are assuming that the ancilla oscillator relaxes faster than the system oscillator as well as that the non-linear 
dynamics are weak compared to the damping of the ancilla oscillator. For the consistency of the following treatment it 
will also be necessary to have XqiNi/k ~ e^. This requirement follows from the second term (the non- linear coupling 
term) on the right hand side of the master equation. This constraint can be achieved consistent with Eq. H15|l . for 
example, by leaving a finite and choosing Ni,Xqi/k ~ e, a regime of low temperature and moderate non-linearity. 
The approximations are also valid at arbitrary temperature in the limit of strong driving and weak non-linearity such 
that Xqi/k ~ and a ~ hold'^^. Here the scaling of the driving strength is chosen to preserve the measurement 
sensitivity which will scale with Agia/K. In this regime the frequency shift of the system oscillator becomes large. 

As mentioned above, in this displaced frame the state of the readout oscillator is close to a thermal state and we 
expand p in the form 



p ^ Po <^ PNi + Pi <^ a\pNi + p\<^ PN^ai + p2<» a\pNiai 



■ P2' ® afpATi + p\, 



PNi aj 



O 



(16) 



Here the pi,i = 0,1,2... act on the system oscillator and the subscripts indicate orders of magnitude in e. The 
scalings of the different parameters with e have been chosen to guarantee the consistency of the expansion. The 
quantity pN^ is the thermal density matrix for the ancilla, which in terms of the average excitation number A^i is 



PN, = 



E 

n=0 



1 



Ni + 1 \Ni + l 



\n){n\ 



(17) 



This is the steady state of the master equation for an oscillator coupled to a thermal bath with temperature given by 
A'^i We have restricted Eq. H16|l to normal ordered terms using the following identities which can be proved from this 
expression for pni- 



PN^ai 



aipNi = 



Ni + 1 



a{pNi 



PNi ai . 



(18) 
(19) 



These normal ordering identities are the key to generalizing the arguments of Wiseman and Milburn to the finite 
temperature case. Using Tri(pAfi) = 1 and Tri{a[pN^ai) = A^i + 1, it can be seen that the system density matrix 
after tracing out the ancilla state is 



Ps = Tri {p} Po + (Ni + l)p2. 



(20) 



Now we substitute Eq. H16|l into Eq. (|14|l and, using the normal ordering identities, derive equations of motion for 
the operators pi, retaining terms in the evolution of pq and p2 up to second order in e: 



Po = - iXoi 
pi = - i 

P2= - iAoiaJao 



* t t t 

a a^a^pi — ap[al^aQ 



+ 2n{Ni + l)p2 + Cnpo + liO (e^) , 



iAoiaofloapo + iXoia— — —po%ao - npi + kO [e ) , 
iVi + i 



ap{ + a 



2kp2 — iXai 



Ni + 1 



ajao, Po 



+ iX 



01 



^1 t 



(21) 
(22) 

(23) 
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where £oPo refers to the damping of the system oscillator described by the last two terms in Eq. H14() . When k is large, 
the equation for pi is strongly damped and quickly decays to the steady state. So we perform adiabatic elimination 
by setting pi = and obtaining an expression for pi 



Pi 



aaQGoPo 



Ni + 1 



Po%ao 



0\ 



(24) 



Substituting Eq. 124|l into Eqs. (|21l I23|l and using the definition of the reduced density matrix Eq. (|20|1 we find, up 
to second order in e, the master equation for the reduced density matrix 



Ps 



\l, \a\' (27Vi + 1) 



ajao. 



a\a(3, ps 



i|wo + Aoi(|a|' + iVi)} 



a\an, Ps 



V {No + 1) (2aQPsal ~ aJaoPs - Psajao) + i^Nq (2alpsao - PsUoal - aoajps 



(25) 



This is the main result of this section. Note that the effect of the adiabatic elimination has essentially been to replace 
fli by Aoi|a|aJao/K, an indication that by measuring the ancilla oscillations it will be possible to obtain information 
about the system phonon number. 



IV. QUANTUM TRAJECTORIES 

Equation H25(l describes the statistical behavior of the system due to the coupling to the thermal bath and the 
indirect coupling to the ancilla thermal bath and the measurement bath, but does not tell us how the measured current 
reflects the system state, or about the correlations of the system dynamics with particular measurement outcomes. 
In this section we derive an equation of motion for the state of the system conditioned on a particular sequence of 
measurement outcomes. This equation is termed a quantum trajectory^'iS'a'^, and results from continually projecting 
onto eigenstates of the current. Since the current effectively measures phonon number, this measurement process will 
tend to force the system towards a pure number state that is consistent with the measurement current. The time 
scale for this to occur will depend on the coupling of the system to the measurement apparatus, which is in turn 
connected to the sensitivity of the measurement. On the other hand, the coupling of the system to a thermal bath 
will lead it to absorb and emit energy from the bath. Thus, in order to determine which number state the system is 
in and track its evolution, it must be possible to distinguish between one number state and the next in a time that is 
short compared to the time scale over which phonons are absorbed from and emitted into the thermal bath. 

A quantum trajectory is constructed as follows. Over each infinitesimal time interval, the system and the measure- 
ment bath states become weakly entangled via the interaction Hamiltonian. As a consequence, at each time instant, 
the state of the system influences the distribution of the possible values of the current / that may be obtained in the 
measurement. In turn, von Neumann projections of the entangled states allow us to calculate the effect of the mea- 
surement of the current on the system state. The appropriate projection is onto the current eigenvector corresponding 
to the measured current value. This results in a stochastic master equation for the state evolution. To implement the 
quantum trajectory approach we perform a simulation by picking the measurements I{t) from the correct probability 
distribution and following the corresponding evolution of the system state. The I(t) curve produced by such a simu- 
lation is representative of a single experimental run, and is a useful predictor of what the experimentalist might see. 
There will be a signal contribution that reflects the system state, as well as a white noise background arising from 
both thermal and quantum noise. 



A. Description of the measurement 



While quantum trajectories are discussed in general at zero temperature in the quantum optics literature, Wiseman 
has discussed the quantum trajectory equations for homodyne detection at finite temperature (Ref.^^, § 4.4.1). The 
demodulated current that reflects changes in the phase of the ancilla oscillation in our setup is mathematically 
analogous to homodyne detection at finite temperature, and so we can adopt these results here. 

The measurement bath is described by the Boson operators bd^n- Since the measurement bath is assumed to be 
large, the finely spaced modes with a smooth density of states leads to a short memory time, a result known as the 
Markov limit. To exploit this is is useful to introduce a global bath operator which captures the combination of bath 
modes that interact with the ancilla oscillating at frequency cui at time t (see Appendix ^ for the derivation of these 
results) 

B, = \ , : E 5d K) &d.„e-^(— ^)*, (26) 

V27rpd(wi)5'd (^i) „ 
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and has time-local commutation rules in the Markov approximation 

Bt,Bl] =6{t-t'). 



(27) 



The operator Bt should be considered to be a linear combination of Schrodinger picture operators, with the phase 
factors of the coefficients depending on the parameter iSi. In quantum optics this is termed the input field operator 
and, roughly speaking, describes the combination of bath modes interacting with the system at time t. In terms of 
these operators the current Eq. (appropriately scaled to remove proportionality constants) is 



I{t) ^Bt + Bl 



(28) 



and the interaction Hamiltonian between the ancilla and the measurement bath in the interaction picture is 

Hl,{t)^-ih^^{Bta\-Bla,). (29) 

with /i ~ TTQd (wi) l^d {^)\^ the ancilla damping rate coming from the measurement bath coupling as before. 

The idea of the calculation is to consider the interaction of the ancilla with the bath at time t, represented by the 
operator Bt, over a small time interval At. It is supposed that each "element" in the time sequence of the bath Bt 
is initially described by a thermal state. Over the interval At the ancilla and bath states become weakly entangled. 
Measurement of the current {i.e., the bath operator Bt + B\) then finds a value of the current equal to an eigenvalue 
/ of the current operator, with the corresponding eigenstate |/), with a probability distribution P{I) given by the 
density matrix of the entangled state in the usual way. 



P{I)^{I\p{t + At)\I). 

The measurement also projects the density matrix onto the eigenstate |/) 

|/) {I\p{t + M)\I) {I\ 
{I\p{t + M)\I) ■ 



(30) 



(31) 



Since the value of the current measured is a stochastic variable, this projection adds a stochastic component to the 
evolution of the density matrix. 

To follow the evolution over a time At it is useful to introduce the normalized operator 



AS = 



At 



Btdt 



which satisfies the commutation rule 



[Ai?(i),Ai?t (t)] 



(32) 



(33) 



At time t the density matrix representing the ancilla and the segment of the measurement bath represented by A_B (t) 
can be written as a direct product of the system plus ancilla p [t) and bath pb {t) density matrices 



p{t)=p{t)®p^ (<), 

and pb {t) is a thermal state. To lowest order the evolution under the interaction Eq. (|29|l gives 
p{t + dt) = p (t) (g) pb (t) + ys^^VAi Iab^gi - alAB, p (t) ® pb (t)] + O {At) . 



(34) 



(35) 

The second term on the right hand side of this equation is the leading order term in the weak entangling of the state, 
and will lead, after projection, to the stochastic part of the density matrix evolution. Using the AB notation has 

made the O (^VKtj size of this term explicit. To derive the deterministic part of the evolution equation we would 

need to keep the O {At) terms, but since these are already known (the master equation in Lindblad form) we will not 
do this here. 

The scheme is now to project this density matrix onto an eigenstate of AB + AB'' chosen with a probability given 
by p{t + dt). Because of the weak coupling of the bath with the system, this will give a small additional contribution 
(actually proportional to %/ At) to the system density matrix depending on the value of the current measured. Since 
the combination of operators AB + AB'' is just the displacement X of the harmonic oscillator represented by the 
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operator Ai3, this projection is most easily done by first transforming the state of the bath into a Wigner function 
representation (see for exampleSS). At time t, the bath oscillator described by AB{t) is in a thermal state and 
the distribution of X is a Gaussian centered a,t X = and with width 2Ni + 1 = coth (fta;i/2fcBT)- Following 
the evolution of the state shows that at time t + At and to O (^\/Atj the distribution of X after the evolution 
corresponding to the operation Eqs. H31(l and H35() remains Gaussian and with the same width, but now centered 
around •\/2/lTrp(() |ai + a|| At. This means that the variable %/ AtX is a Gaussian random variable given by 



AtX = y^lai + a\ )At+ y/2Ni + IdW, 



(36) 

with dW a Wiener increment with dW^ — At. These results give us expressions for the measured current and the 
effect of the measurement on the system density matrix. 

Since the current is XjyJ At^ the first important result is that the measured current integrated over time At is 



/ (<) At = (ai+ a\ )At + y/2Ni + IdW, 



or in differential form 



/ it) ^^l^ai+ a\) (t) + V2iVi + 1^ (t) 
where f (t) — dW/dt represents white noise with correlations 

(eW> = o, 

{at)Ht')) = Ht-t'). 



(37) 



(38) 



(39) 
(40) 



The second result is for the increment of the system density matrix after evolution through At and projection by the 
measurement (cf.^^ Eq. (4.113)), 

dp^' {t) = {X\p{t + At)\X)/p{X) 

{Ni + 1) (aip^* + p^'flt) - iVi (atp^' + p'^^al) -Tr i^aip^' + p''a\} p''*] + O (At) . (41) 



AtX- 



2p, 



2Ni + 1 



Replacing the stochastic variable X by the expression Eq. 1^^. and retaining only the O ^VAt^ term gives 

(A^i + l)(ai/9''' + p'^'al) - Ni{a\p'' + p'^'ai) - (ai + 4)^^*1 dW + O (At) . 



dp"" 



2p 



l + 2Ni 



(42) 



Equation (|42|l is the stochastic term that must be added to the density matrix evolution of equation (|13|1 to give us 
the stochastic master equation for the density matrix conditioned on the measurement outcome. Note that the noise 
term dW appearing in Eq. H42|) is the same as that appearing in Eq. H37|l . so that it is related to the actual current 
measured I{t) by Eq. 



dW = (l{t) - y%(ai+al)) Ai/Vl + 2iVi. 



(43) 



B. Adiabatic elimination on the stochastic master equation 

Just as we did for the master equation it is possible to adiabatically eliminate the ancilla coordinates and find a 
stochastic master equation for the system alone. Using the same expansion for the conditioned density matrix Eq. 
(|16|l we can determine stochastic equations for pf from Eq. (|42|l . We obtain the set of differential equations 

dpf = /^^^ {(^1 + 1) {pf + Pf^) {pf + pf) Po} dW + 0{dt\ (44) 
dpf - /^J^ {(^1 + 1) (2P2^ + pf) - {pf + pf) pf] dW, (45) 
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We have written only the stochastic contributions; the terms proportional to dt are the same as in the adiabatic 
elimination on the ordinary master equation. 

Now we will do the adiabatic elimination as we did for the deterministic master equation. As before we wish to set 
dpf to zero. However, since dpf is stochastically driven this will not be precisely true even in the steady state. In 
order to replace pf by its mean value at long times it is necessary to estimate the size of the fluctuations resulting 
from the dW-teim. Following the analysis of Doherty and Jacobs^S we integrate the stochastic term over the decay 
time of the ancilla oscillator and compare its root mean square magnitude with the deterministic terms. Consider the 
full equation for dpi 



dpf 



iXoia 
npidt + 



t St 

%aQPo 



Ni .t t ■ 



Ni + 1 



dt 



y^^^ {(TVi + 1) {2pf, + p,) {pf + pf) pf} dW + kO {e^) dt. (47) 



We integrate this over a time At ^ l/n and use the fact that the mean values of po,pi must be slowly varying over 
this time to obtain 



At 



dpi 



iXoia 



alaopo 



A^i + l 



At - KpiAt 



+ + 1) (2P2' + P2) - {Pi + Pit) Pi}J^^' dW{t') + O (e3) 



(48) 



(49) 



The random number AW = /g^* dW{t') is Gaussian distributed with mean zero and variance At ~ l/^*^, thus the 
root mean square size of AW is As a result the stochastic term in the update of pi scales like e^^^ and is 

negligible in comparison to the deterministic terms which scale like e. As a result Eq. 1)24(1 holds exactly as before. 
Using Eq. (|20|l and ia — ~ia* = \a\ we finally obtain the stochastic master equation (SME) for the system 



dps = - 



Xli \a\' (27Vi + 1) 



ajao, Ps 



dt 



^{cJo + Aol(|a|' + iVl)} 



ajflo, Ps 



dt 



V (Nq + 1) (2aoPsal ~ aJaoPs - Psfljao^ dt + vNq {2alp^ao - p^aoal - a^alp^ dt 



alaoPs + PsCilao - 2{alao)ps 



dW, 



where 



(50) 



k = fiXl^lal"^ / {2Ni + I) k'^ 



(51) 



Again the noise dW is related to the measured current, which using Eq. H24|l can be written in terms of the system 
phonon number 



I{t)At = yj2Ni + 1 {2\/2k (alao^ At + dW 



(52) 



V. RESULTS 



To further understand the consequences of the stochastic density matrix equation Eq. H5U|) we consider a case in 
which the initial state is a mixture of number states, ps = '^nPn\n){n\. (A thermal state is an example of such a 
state.) The solution of the stochastic master equation, Eq. (I50() . remains a mixture of number states if the initial 
state is a mixture of number states. For such an initial state the stochastic master equation can be reduced to an 
equation for the weights p„, which takes the form 

dpn = -2v{No + l)[npn - {n + l)pn+i]dt 

- 2iyNo[in + l)p„ - npn-i]dt - 2^/2fc - '^'^"'^ P^'^^' 
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FIG. 2: Plot of the solution of Eq. 15311 without the stochastic component, A; = 0, with the initial state |1) (solid line) and |2) 
(dashed line). 



Since mixtures of number states are invariant under changes of phase and the number states are eigenstates of the 
Hamiltonian, neither the phase diffusion term nor the Hamiltonian terms in the stochastic master equation contribute 
to the evolution of the phonon number distribution. It can also be shown that this system of equations also describes 
the evolution of the phonon number distribution p„ = (n|p|n) for an arbitrary (not necessarily diagonal) initial state 
P(0). 

Equation (|53|l is our central result for analyzing the behavior of the measurement protocol. The first two terms of 
Eq. H53I) containing dt describe emission into and absorption from the thermal bath coupled to the system. We will 
see that the second, stochastic, term tends to concentrate the distribution p{n) onto a single value of n. The two 
effects are characterized by the two time scales, and k~^, and the resulting behavior depends on the ratio of these 

two times. We will discuss the competition by calculating the occupation number of the system, ^ajao^ (t), which is 

given in terms of p„ from numerical simulations of Eq. (|53|l by 



(<) = ^np„ it). 



(54) 



Firstly we turn off the stochastic component and consider the solutions given by the deterministic part of Eq. 
(|53|l . Figure 121 is a plot for fc = starting from two different initial states |1) and |2) and for a bath temperature 
corresponding to an average occupation number A^o = 1.62. The plot shows that the deterministic terms in Eq. H53() 

drive the system towards a mixed (thermal) state, so that the ensemble average of (^ajao^ (t) gradually reaches the 

thermal average at the bath temperature. This is true regardless of the initial state. Note that the deterministic part 
of Eq. H53|) also describes the average over all measurement outcomes even for nonzero fc, since the stochastic term 
averages to zero. We can define the characteristic time the system resides in a given number state, which we call the 
dwell time tdwcii, as the reciprocal of the initial transition rate given by Eq. (|53|l with fc = 



^dwcll 



dpn 

dt 



1 



P„(0) = 1 



2iy[No{n + l) + {No + l)n] 



(55) 



Note that the dwell time depends on the initial state n, and also on the temperature of the bath through Nq. 

We turn now to the dynamics resulting from the measurement process in the absence of coupling to the thermal 

bath. Figure 13 shows results for ^ajap^ (t) for a simulation of Eq. H53|l with v ~ and an initial condition of a 

thermal state. Figure^ shows the individual probabilities pn for n = 0,1,2,3 for the same simulation. All number 
states are present initially, but eventually the system is projected onto state |1) in this simulation. In other runs, 
with different random numbers for the stochastic term, different final states result, as expected. The plots show that 
the stochastic term tends to project the system state onto a pure number state on a time scale of order k~^. We call 
this time the collapse time ^coU- Since no coupling to the thermal bath is present in these simulations, once projected 
onto a number state, the state is stationary. The collapse onto a number state can actually be shown analytically 
using the solution of the system of equations (|53|l due to Jacobs and Knight^. 
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FIG. 4: Plot of Pnit) for a simulation of Eq. I53II with 1^ = for the states |0) , |1) , |2) , |3). The initial state is a thermal state 
with the average occupation number 1.63. The figure is for the same simulation as in Fig. |21 
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FIG. 5: The evolution of the phonon number ^ajaoy (t) given by Eq. H53^ using No — 1.62 and k/iy = 250 (first panel) and 
k/i/ = 5 (second panel). 
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FIG. 6: Histogram of (alao) {t) for a long simulation with k/u — 150 and No = 1.62 (first panel) and k/u = 15 (second panel). 



For the phonon number ^ajaoy {t) to take on discrete values with both the thermalization by the coupling to the 

bath and the projection by the measurement process present, we need tdwcii ^ ^coii- This is illustrated in Figs. |S1 
which show results for the cases tdweii ^ ^coii and the ^dwcii ^ ^coii with the fixed value of A^o = 1-62. We use 
values of k/i' = 250 giving tdweii/icoii = 153 for state |0) and tdwon/^coii = 42.4 for state |1), and k/iy = 5 giving 
idwcii/icoU = 3.06 for state |0) and idwcii/^coii = 0.85 for state |1). The jumps in the occupation number are clearly 
evident in the former case, but are not seen in the latter case. The discreteness in the phonon number is shown more 

clearly by plotting histograms of (ajao) (i), Figs.El again using a fixed value of A'o = 1.62 but with different values of 



k/v equal to 150 and 15. (A bin width A ^ajooy — 0.1 is used.) The clustering of the (^ajaoy (t) values around integral 

values is clearly evident for k/i' = 150, is still identifiable for k/iy — 15, and completely absent for k/f = 3. The 
increasing sharpness of the jumps with larger k/v can be seen in a more quantitative manner by plotting the standard 



deviation of the phonon number from integer values, the time and ensemble average of 

as a function of k/i^ (see Fig. (T)). 

Since tdweii is dependent on the temperature, the condition id 



(t) - Int 



ajao 



it) 



well ^ icon effectively places a limit on the temperature 
of the system oscillator even for large k. Setting n ~ Nq in Eq. for idwcii, this inequality gives the condition on 
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FIG. 7: The deviation of <^aQaoy (t) from integral values as a function of k/u. The deviation is defined as the average of 
ajao^ (t) - Int (^ajfflo^ • 
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FIG. 8: Evolution of (^ajaoy (t) at a temperature corresponding to A'o = 1.62 (first panel) and A^o = 20 (second panel), with 
vNo/k — 0.0108 and from an initial state |2). 



the temperature for jumps in the number to be seen 



lyNoiNo + I) /2k < 1. 



(56) 



In order to keep the same resolution for observing clear jumps as at low temperature, k/iy must be increased as 
temperature increases. This is not an easy task for the experimenters: for example, an oscillator with 1 GHz resonant 
frequency, which is the highest frequency currently reported for a mesoscopic oscillator—, at T = 0.1 K the average 
occupation number is Nq — 1.62. When the temperature is raised to T = 1 K, the value rises to A^o = 20. Thus if we 
demand the same resolution for jumping in both cases, the sensitivity of the measurement at the higher temperature 

must be increased by a large factor. This is illustrated by Figs. |H1 which show ^ajag^ (t) over time for different 

temperatures corresponding to iVo = 1.62 for k/v = 150, and Nq = 20 for k/v = 1850. The product vN^/k has 
been kept constant at 0.0108 in order to provide the same resolution for the jumps. Also notice from Eq. (|55|l that 
^dwcii decreases with the system state n making it difficult to recognize the discrete jumps when the system state is 
at higher n. 

We have so far considered the possibility of observing discrete occupation numbers in terms of the behavior of the 
variable ( ajao ) (t). In actual practice the occupation number must be inferred from the measured current I{t), and 
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FIG. 9: Filtered current using a running average over various window sizes, for parameters k/v = 250 and A^o = 1.62. The dotted 

line is ^ajoo^ {i) given by the stochastic density matrix. The current was first averaged over a time interval of /sAt = 0.15. 

Upper left: current observed; upper right: window size feAt = 4.5; lower left: window size A;At = 7.5; lower right: window size 
kM = 10.5. 



is obscured by the noise in this variable. A simple scheme to reduce the effect of the noise is to average the signal 
over a sliding window. We can define the measurement time tm as the averaging time required to give unit signal to 

noise ratio. Thus we equate the signal S, given by averaging the current for unit phonon number ^ajao^ = 1 over 

the measurement time 



with the noise N over this averaging time 



N : 



Setting S/N = 1 gives the measurement time 




V'2iVi + It 



1 

8k' 



1/2 



(57) 



(58) 



(59) 



For jumps in the measured phonon number to be detected in the current we would need tdwcU ^ ^m- Notice that the 
measurement time and the collapse time are comparable. This means that if the experimenter can infer the system 
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number state through the measurement current, then the system is actually projected to that state on the same time 
scale. The results for different averaging times At is shown in Fig. O For kAt equal to 4.5 or 7.5 the averaging is 

sufficient to display the steps in ^ajap^ (t) without too much rounding of the transitions. 

The simple averaging is not actually the optimal way to extract ^agao^ (t) from I{t). In principle a better approach 

is to use the stochastic master equation to reconstruct the dynamics of the system given the initial ensemble Ps(^o) 
and the measured current I{t). This can be seen more readily if we rewrite Eq. H52|) as 

dW^-j=^=={lit)-2V2k(^alao){t)]dt. (60) 

In our simulations we draw I{t) at random from the appropriate distribution and find the stochastic density matrix 
Ps- However, using I{t) from experimental data, the experimenter can in principle propagate Eq. 1)53(1 using Eq. 1(60(1 
and then can estimate the phonon number at each time from Eq. ((54(1 . This procedure is itself a low pass filtering that 
reduces the noise on the measurement current, corresponding to the optimal filtering for our model of the systemMi^ 



VI. PARAMETERS AND CONSTRAINTS 



The results of the previous section show that a large value of the ratio k/v is crucial. To analyze the interplay of the 
parameters of an experimental realization, we simplify the expression for k/v by assuming that most of the damping 
of the ancilla comes from the necessary coupling to the measurement device, rather than from the extra thermal bath, 
i.e. /X ~ K. Then using the expression of k from Eq. ((51(1 and v =^ luo/2Qo we get 

^^4(27Vi + l)-iQoOi-f— Vlap. (61) 

Equation 1(611) shows that the success of the measurement procedure is favored by large oscillator quality factors, 
large driven response |a|, and a large value of the anharmonicity coupling factor Aqi/cji. In addition, as we have seen, 
detecting individual jumps becomes harder as the temperature increases. Increasing the quality factors of mesoscopic 
oscillators is an active area of research. Currently, values of order 10'^ to 10^ seem possible. If these could be raised to 
the values characteristic of more macroscopic oscillators of the same material, of order 10^ or even higher, the detection 
of individual phonons would become correspondingly easier. The frequency ratio loi/ujo appearing in Eq. 1(61(1 must 
be less than unity for our detection scheme, but will probably not be too small because of geometry constraints. Thus 
the main parameters available to optimize the experimental geometry are the anharmonicity factor Aqi/wi and the 
dimensionless measure of the driven displacement of the ancilla, |a| (the number of phonons in the driven state). We 
now consider these factors in more detail. 



A. Anharmonicity Coefficient 

The interaction Hamiltonian for the system and ancilla oscillators Eq. ((TJ can be written 

^RWA _ fiujQalao + h[uji + Aoi^o] a\ai, (62) 

with no the system phonon number. This equation implies that Aoi can be estimated as the frequency shift of the 
ancilla oscillator for a single quantum (no = 1) of the system oscillator. 

For the prototype geometries using the two orthogonal flexing modes of a single beam, or parallel flexing modes 
of two longitudinally coupled beams, the nonlinear coupling arises from geometrical effects. At second order, the 
transverse displacement in one mode gives a longitudinal strain, which then changes the frequency of the second 
mode. The strain generated by the flexing motion and the frequency shift associated with this strain can be derived 
using elasticity theory, and have been calculated by Harrington and Roukes^. A demonstration of such frequency shift 
detection, and direct measurement of Aoi between two coupled beam has recently been reported^. The longitudinal 
strain produced by a single quantum in the fundamental flexing motion is 
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where toq is the mass and Lq is the length of the system beam. Then the ancilla frequency shift caused by this strain 
is 

Aoi=a;i^X^, (64) 

where C is a geometric factor = 3 for clamped beam boundary conditions) and Li, di are the length and thickness 
of the ancilla beam, respectively. Introducing a dimensionless quantity, 

fi^ 1 

(65) 

rriidf rujji 



then the scaled coupling coefficient can be expressed as 

Aqi __ C miUJiLl 



R. (66) 



Since the factor of the ratio of the two mode parameters will not be too large or small, the most important quantity 
determining the anharmonicity factor, which must not be too small for the success of the measurement scheme, is the 
dimensionless ratio R. This will typically be a small number. The need for small devices is seen from the scaling of 
this parameter with the dimensions. 



B. Driving strength 

The detection scheme we have considered is to measure the phase of the driven response of the ancilla oscillator. 
Since the detection scheme is magnetomotive, it is natural to consider the use of magnetic driving in estimating 
the size of the displacement parameter \a\. For magnetic driving using a current /drive in a magnetic field B the 
dimensionless displacement can be estimated as 

H = Ql^%^^/i^. (67) 

Again the important role of R in limiting the size of |a| in this analysis is apparent. 

We must also recognize that the size of |a| might be limited by other physical constraints, rather than by the 
available drive strength. One constraint might be to avoid undesired nonlinear effects in the driven beam itself. For a 
classical oscillator, at sufficiently large drive amplitudes nonlinear frequency pulling leads to a multiplicity of solutions 
and instability. This occurs when the nonlinear frequency shift is comparable with the width of the resonance coi/Qi. 
Using the same type of estimate for the nonlinear frequency shift as in Eq. (|64|l shows that this occurs for 



> 



VOIR 

A more detailed, quantum mechanical analysis of the driven nonlinear oscillator will be presented elsewhere^^. 



(68) 



C. Example Configuration 

As a first estimate of the order of magnitude of the quantities introduced above we will construct an example 
configuration using parameters that seem plausible with current technology. 

Recently, oscillators with resonant frequencies as high as 1 GHz have been fabricatec^ using silicon carbide. Thus we 
consider two flexing modes with resonant frequencies of too — 2.3 GHz and uji — 0.36 GHz so that ujq — toi ^ Aqi, i', k 
are satisfied. For this value of the system oscillator frequency, hw^/k^T is unity at a temperature of about 0.1 K. The 
oscillators in Ref.^ were not very small, but it is expected that the structure can be scaled down while maintaining 
the high oscillation frequency. We therefore suppose smaller dimensions consistent with these frequencies, namely 
dimensions are 0.6 /im x 0.04 /im x 0.07 /zm for the system beam and 0.6 /Ltni x 0.04 /itm x 0.01 /Ltni for the ancilla 
beam. With these parameters we obtain 



i? = 4.26x 10~^ 



(69) 
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The factor R occurs squared in k/v (via Eq. H66|l . which is required to be large, and so this small factor must be 
mitigated by the other quantities in Eq. H66|) . i.e. large values of Q and a large driven amplitude \a\. Suppose the 
Q of the system oscillator is Qo = 10000 and Qi = 1000. For the size of |a| first consider the magnetic driving. A 
magnetic field B = 10 Tesla and /drive — 1 /^A can raise the driven response to \a\ ~ 10^. To reach k/i/ ^ 1 the 
nonlinear coupling Aqi/wi that is required is then 

Aoi/cJi = 4.9 X 10"^ (70) 

With the given beam dimensions and the geometric nonlinearity, the anharmonic coupling coefficient is actually 
Aoi/wi — 1.3 X 10^^^, about three orders of magnitude smaller. One possible way to increase this value might be 
to engineer the geometry of the oscillator so that the anharmonic coupling is larger than in the simple geometric 
nonlinearity we have considered'^^ . Another way to increase k/iy is to increase |q;| using a different driving scheme, 
although for the value or R in Eq. 1)69(1 the limit in Eq. ((68|l is already exceeded for \a\ > 10'^ , so that engineering 
the geometry to reduce the self- nonlinearity might be necessary. An obvious way to increase k/v to values greater 
than unity is to use oscillators with smaller dimensions, for example carbon nanotubes. 

VII. CONCLUSION 

We have analyzed a scheme to observe quantum transitions of a mesoscopic mechanical oscillator. The non-linear 
coupling shifts the frequency of a second (ancilla) oscillator proportionally to the excitation level of the first (system) 
oscillator. This frequency shift may be detected as a phase shift of the ancilla oscillation when driven on resonance. 
In principle, a QND measurement is possible if the coupling constant between the two oscillators Aqi is much smaller 
than the resonance frequencies of the oscillators, as will usually be the case. We have derived the master equation for 
the system density matrix first integrating out the environment and measurement degrees of freedom, and then by 
removing the ancilla operator using the fact that the time scale of the system and ancilla dynamics are quite different. 
The master equation has three components: phase diffusion as a result of the measurement backaction; a constant 
energy shift due to the excitation of the ancilla oscillator, and number state transitions due to the interaction with 
the thermal bath (the environment). 

The measurement process introduces a stochastic component into the system dynamics and we have obtained 
the stochastic master equation corresponding to our measurement scheme. From the stochastic master equation we 
identify two competing tendencies that can be characterized by two parameters. One is the coupling strength i' 
of the system and thermal bath, which is associated with the dwell time tdweii between transitions. The other is 
the coefhcient fc, associated with measurements, which includes not only the coupling strength of the system to the 
measurement bath but also the anharmonic coupling strength between the oscillators, the driving amplitude. This 
coefficient is related to the measurement time t^^ that is needed for a measurement to be able to produce an outcome 
with certainty. To observe clear quantum jumps we would need idwoU S> tm- If this condition is not satisfied, then 
the experimenter cannot infer the energy eigenstate of the system from the observed current. 

Although our simple estimates based on plausible lithographically prepared oscillators yield values for the ratio 
tdwcii/tm too small for the observation of individual phonons, enhancements to the geometry and the trend to smaller 
device sizes should improve the outlook. The basic scheme and theoretical techniques developed here are fairly general, 
and in particular are not restricted to zero temperature, and so can be also used for other applications such as single 
spin detection and noise analysis for a solid state based quantum computer. Such possibilities might open up a new 
stage for observing quantum dynamics in mesoscopic systems. 
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APPENDIX A: BATH FIELD OPERATORS 



In this appendix we describe more fully the time-local measurement bath operators Bt introduced in section FlVAl 
The description in terms of finely spaced modes of the bath with a smooth density of states leads to the short memory 
or Markov property of the bath, that can be expressed in terms of the time-local commutation rules for Bf. In the 
main text we introduced the global bath operator as Eq. 



Bf 



{uJn —LJl}t 



^/2^^pd{uJl)gd (uji) 

We first derive the commutation rule Eq. H27|l. Substituting Eq. (|A1|) in the commutator gives 



[Bit),B^ it')] 



27rpd (t^i) [gd (wi)]' 



(Al) 



(A2) 



and using 



[B {t) , fit (t')] ^ _ 

27rpd (^i) [gd (wi)] 



Changing the sum to integral form 



E 



duj pd (uj) 



gives 



[B{t),B^ (0] - 



27rpd (wi) [gd (wi)] "'0 



(A3) 



(A4) 



(A5) 



Since pd (w) and gd{^) are slowly varying functions around the ancilla oscillation frequency lj = uji we can approximate 
these as pd (w) ~ pd (wi) , gd{^) — gd{^i) and pull them outside of the integral. Then introducing e = uj — uji and 
extending the lower range of the integration over e to ~oo leads to the desired result 



[B (t) , Bt {t')] - ^ (*-*') =5{t-t'). 

The interaction Hamiltonian for the ancilla oscillator and the measurement bath is, from Eq. H12|l . 



(A6) 



(A7) 



where we have moved to the interaction picture with ai{t) = aie^*'^^* and bd,n{t) — &d,rie^*'^"* the ancilla and bath 
operators in this picture, and gd (^n) is the coupling strength. The interaction Hamiltonian can be written in terms 
of the bath operators Bt as 



where the coefficient p is 



p = TTQd (wi) \gd {i^i)Y 



(A8) 



(A9) 



as before, and we have used the fact that the ancilla interacts predominantly with bath modes near frequency oji 
and again have assumed a smooth variation of the density of states and coupling constant so that we can make the 
replacements pd {uj) ~ pd (wi), and gd{uj) — gdiuJi)- 
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